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ABSTRACT 

The association of an electromagnetic signal with the merger of a pair of supermassive black holes 
would have many important implications. For example, it would provide new information about gas 
and magnetic field interactions in dynamical spacetimes as well as a combination of redshift and 
luminosity dist ance that would enable precise cosmological tests. A proposal first made by |Bode "&~ 



Phinney (20071) is that because radiation of gravitational waves during the final inspiral and merger of 



the holes is abrupt and decreases the mass of the central object by a few percent, there will be waves in 
the disk that can steepen into shocks and thus increase the disk luminosity in a characteristic way. We 
evaluate this process analytically and numerically. We find that shocks only occur when the fractional 
mass loss exceeds the half-thickness (h/r) of the disk, hence significant energy release only occurs for 
geometrically thin disks which are thus at low Eddington ratios. This strongly limits the effective 
energy release, and in fact our simulations show that the natural variations in disk luminosity are 
likely to obscure this effect entirely. However, we demonstrate that the reduction of luminosity caused 
by the retreat of the inner edge of the disk following mass loss is potentially detectable. This decrease 
occurs even if the disk is geometrically thick, and lasts for a duration on the order of the viscous 
time of the modified disk. Observationally, the best prospect for detection would be a sensitive future 
X-ray instrument with a field of view of on the order of a square degree, or possibly a wide-field radio 
array such as the Square Kilometer Array, if the disk changes produce or interrupt radio emission 
from a jet. 

Subject headings: accretion, accretion disks - black hole physics - galaxies: nuclei - gravitational 
waves — relativity 



1. INTRODUCTION 

There has recently been significant interest in whether 
the coalescence of two supermassive black holes can pro- 
duce an identifiable electromagnetic (EM) signature. If 
a redshift could be extracted from such a signature (via 
localization of the host galaxy), then it would provide 
a powerful cosmological probe when combined with the 
expected estimate of the luminosity distance, which will 
have a precision of ~ 10% mainly limi ted by weak gravi- 
tational lensing (Hughes & Holz|2003). It would also al- 



in the 



low precise tests of whether gravitons travel at the speed 
of light, as required by general relativity. 

While the merger itself produces no electromagnetic 
emission, if there are significant electromagnetic fields or 
mass nearby in an accretion disk then there are various 
possibilities. Milosavljevic & Phinney (2005) point out 
that for some disk accretion rates ana binary mass ra- 
tios, the binary reaches a point in its coalescence such 
that further inspiral by emission of gravitational waves 
occurs more rapidly than the disk diffuses inwards. This 
leads to a hole in the disk which is filled gradually af- 
ter merger, leading to a source that brightens over weeks 
to years dep ending on various parameters. Several re- 
cent authors ([Shields & Bonning | 2008 [ |Lippai et al.| 2008 



Schnittman & Krolik 20081) discuss consequences of the 



recoils from asymmetric emission of gravitational waves 
during the coalescence, from prompt shocks to delayed 
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emission lasting mil lions of years. Emission might occur 
l ate inspiral ([A rmitage & Nat arajanp 
b||2008| |Haiman et al.||2008|) from effi 
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fc Loeb 2008 Haiman et al. 2008) from effects such as 
enhanced accretion, periodic Newtonian perturbations, 
or shearing of the disk due to g ravitational waves. Ear- 



lier precursors are also possible (|Dotti et al.|2006; Kocsis 



et al. 2008 Haiman et al. 2009), and in some cases the 
error volume from the gravitational wave signal may be 
small enough that the host galaxy can be identified by 
its morphology or mass, or by the presence of an Active 
Galactic Nucleus (AGN) (e.g., |Kocsis et al. 120061 ). 

A suggestion first made by |Bode fc Phinney (2007) 
is that the effectively instantaneous mass-energy loss to 
gravitational radiation would leave the disk in a non- 
equilibrium state, and that adjustment of the disk could 
produce observable radiation. Their specific suggestion 
was that because fluid orbits would be elliptical after 
mass loss, circularization would release extra energy that 
could be detected. In this paper, we examine this sce- 
nario and find that in most circumstances this energy 
release is too small, and the timescale too large, for such 
emission to be realistically detectable. We show that in 
principle, for systems accreting at a high enough rate 
that the inner edge can stay within a factor of ~ 2 of 
the binary semimajor axis all the way down to binary 
dynamical instability (i.e., when the black holes begin 
to plunge into one another), relativistic effects enhance 
the energy release significantly. Rapid release of this en- 
ergy, however, requires the development of shocks. As 
we demonstrate, shock formation only occurs if the disk 
is sufficiently geometrically thin. Such disks are char- 
acterized by a long viscous timescale and thus are not 
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sufficiently relativistic at the point of merger for there to 
be substantial energy release due to circularization. 

Here we show that a stronger signal may be obtained 
from the prompt drop in emission that occurs because 
the inner edge of the disk retreats after mass loss and 
then takes some time to fill in. The resulting flux deficit 
compared to the emission prior to mass loss is in the de- 
tectable range for current instruments such as Chandra, 
and observable for planned observatories such as the In- 
ternational X-ray Observatory, if the angular localization 
of the gravitational wave event prior to merger is suffi- 
ciently precise. In § 2 we discuss the basic energetics and 
timescales for both Newtonian and relativistic disks. In 
§ 3 we give the results of our hydrodynamic and mag- 
netohydrodynamic simulations. In § 4 we discuss the 
potential observability of this effect, particularly when 
measured against the backdrop of galaxy mergers and 
AGN variability, and in § 5 we give our conclusions. 

2. BASIC ENERGETICS FOR NEWTONIAN AND 
RELATIVISTIC DISKS 

2.1. Fractional mass loss 

Our starting point is the typical fraction of total binary 
mass lost to gravitat i onal r adiation during a coalescence. 
|Tichy fc Marrone tti ( 2008 ) analyze the results of numer- 
ous simulations and propose an empirical relation that 
accurately estimates Ara/M tot , where Am is the mass- 
equivalent lost to radiation and M tot was the mass of the 
binary when the black holes were well separated. They 
find 



Am/Mtot « 0.0485(4/?) + 0.013(a, + b z )(4 V ) 2 



(1) 



Here, for a black hole binary of masses m a and ra&, rj is 
the symmetric mass ratio n = m a m5/(m a + m^) 2 . This 
has a maximum of n max = 1/4 for an equal-mass bi- 
nary. We also define a z and b z to be the dimensionless 
angular momenta of the two holes in the direction of 
the orbital angular momentum, a z = cJ ajZ /(Gm 2 l ) and 
b z = cJ byZ /(Gml). 

For gas-rich mergers such as those li kely to produce the 
strongest electromagnetic signatures, Bogdano vic et al. 
(2007) argue that torques from the gas would tend to 
align the spins of the individual holes with the orbital 
axis and hence with each other. Accretion is thus usu- 
ally prograde, hence black holes in gas-rich environments 
accreting for extended times are expected to spin rapidly; 
observational evidence for high spi ns has been presented 
via analysis of iron Ka lines (e.g., Iwasawa et al. 1996 
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& Reynolds 2006). For these objects a z « b z w 1 and 
the mass loss is maximized. A reasonable approxima- 
tion to the fractional mass loss is therefore Am/M tot « 
0.05(4^) + 0.025(4t7) 2 . This is roughly 7.5% for equal- 
mass systems, 5% for a 1:3 mass ratio, and 2% for a 1:10 
mass ratio. 

A useful simplification is that most of the mass is lost 
at the very end of coalescence, hence from the perspective 
of the disk there is a rapid reduction of the mass of the 
central object. For example, when two equal-mass non- 
spinning black holes coalesce, roughly 70% of the total 
energy radiated com es out in the fin al orbit plus merger 
and ringdown (e.g., |Pretorius||2005| and numerous sub- 
sequent papers). We will therefore treat all such losses 
as instantaneous. In addition, we consider the effects of 



mass loss to be decoupled from the gravitational wave 
kick that a remnant black hole may receive after coa- 
lescence. This is justified since the orbital speed in the 
innermost part of the disk is large compared to the ve- 
locity of the kicked black hole. Consequently, any EM 
signatures associated with the kick occur primarily at 
much larger radii, and hence timescales, than the inner 
regions of interest here. 

2.2. Newtonian energetics and timescales 

For simplicity, we will consider first a disk with an 
inner edge far enough from the merging binary that a 
Newtonian treatment around a point mass is justified. 
Suppose that the fractional mass loss is e <C 1, i.e., the 
mass after merger is M = Mq(1 — e) where Mq is the 
initial binary mass. We also presume that initially the 
disk gas moves in essentially circular orbits, meaning that 
the inward radial speed is much smaller than the orbital 
speed. Our final assumption is that the disk is initially 
axisymmetric as opposed to having significant azimuthal 
variations due to spiral density waves induced by the 
binary. 

A fluid element at radius r then has initial specific 
angular momentum £q = ^/GM^r. After reduction of 
the central mass, the specific angular momentum is un- 
changed but the fluid element is now at the pericenter 
of its orbit. If we now consider fluid elements all at the 
same radius as constituting an annulus, then the annu- 
lus will oscillate in phase. Other annuli will do the same 
thing but at different frequencies, meaning that to lowest 
order the accelerations on an annulus due to other annuli 
will all be radial and thus not change its specific angular 
momentum. We therefore assume that the specific angu- 
lar momentum in a given fluid element is a constant and 
then minimize energy under this constraint. 

The initial specific energy after mass loss is E tot = 
\v 2 - GM/r = GM /(2r) - GM (1 - e)/r. Rewriting, 
Etot — — GMq(1 — 2e)/(2r). The initial semimajor axis is 
then given by E tot = -GM/(2a) = -GM (1 - e)/(2a). 
Equating these two expressions for E tot and keeping only 
the lowest-order terms, r « (1 — e)a. Thus, the eccentric- 
ity is e « e. The specific angular momentum of an orbit 
with eccentricity e is /(e) = y/GMa(l — e 2 ). The change 
in binding energy as a result of orbital circularization is 
given by AE = [/(e) 2 - /(0) 2 ]/(2a 2 ) = -GMe 2 /(2a), 
so the fractional change in the binding energy becomes 
AE/E tot = e 2 = e 2 . 

The timescale on which this energy is released is un- 
likely to be shorter than the time needed for neighboring 
annuli to develop significant relative radial motion. Or- 
bits with semimajor axes differing by a factor ~ (1 + e) 
will intersect each other when their radial epicyclic oscil- 
lations get out of phase by ~ 1 radian. Since in Newto- 
nian gravity the radial epicyclic frequency n equals the 
orbital frequency, the minimum timescale at radius r be- 
comes 
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where P O rt>(0 is the orbital period at r. The actual time 
of energy release could be significantly greater than T m i n 
if, for example, no shocks ever develop or there is a delay 
in radiating the energy from the disk due to large vertical 
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optical depths ( Schnittman fc Krolik|20 08). This in turn 
implies a maximum specific luminosity of 



AE/T n 



10e 3 E tot /P olh (r) . 



(3) 



This is to be compared with the natural energy release 
of the disk as it flows inwards. For example, consider 
a fairly large mass loss of e = 0.05. Th en L TT1Ry ^ 
10~ 3 E tot / P or b (r) . For a Shakura-Sunyaev (jShakura 



Sunyaev] 1973) viscosity parameter a = 0.1 and disk 
half-thickness h/r = 0.1, classical disk theory predicts 
an inward radial speed of v r = |tt(^/r) 2 v or b, where 
v rb is the orbital speed. The characteristic time to 
move from a radius r to a radius r/2 is thus roughly 
^orb( r )/[67ra(/i/r) 2 ] w 100P or b(^), so the specific lumi- 
nosity is L ~ 10 _2 ^tot/^orb(^)- The energy release from 
mass loss would in this case be only a ~ 10% perturba- 
tion. Physically, of course, it is the value of h/r that is 
determined by the luminosity of the disk and not vice- 
versa, but it is still convenient to parametrize the lumi- 
nosity in terms of the aspect ratio. While this estimate 
of a ~ 10% perturbation applies to a steady disk with 
the given parameters, we further note that disks with 
lower values of a and h/r would feature even smaller 
perturbations. 

2.3. Relativistic energetics 

If the inward radial speed of the disk is sufficient to 
keep up with the inspiraling binary, then near merger the 
disk experiences relativistic dynamics. As we discuss be- 
low, it is unlikely that in such a case shocks would be pro- 
duced after the mass loss, because rapid inward motion 
requires a thick disk with h/r ~ 1, hence the disk sound 
speed, c s ~ (ft/r)^, is significantly higher than the ra- 
dial oscillation speed, v r ~ ev^, after merger. Nonethe- 
less, for completeness we show that if somehow angular 
momentum transport is strongly enhanced (thereby al- 
lowing a thin disk to follow the binary all the way to 
merger), the resulting energy release could be boosted 
by large factors. 

A key feature of relativistic dynamics is that the spe- 
cific angular momentum of circular orbits reaches a mini- 
mum, at the innermost stable circular orbit (ISCO). The 
specific angular momentum curve is thus nearly flat near 
the ISCO, hence the sudden loss of mass at merger can 
lead to much greater motion of the accreting gas than in 
the Newtonian case. This in turn produces larger and 
potentially faster energy release, thus possibly greater 
observability. 

The spacetime will not be Schwarzschild, both because 
the source of gravity is actually a binary and because 
rapid rotation introduces frame-dragging. Nonetheless, 
the ISCO is a general feature and therefore we present a 
calculation in the Schwarzschild spacetime as represen- 
tative of the expected effects. 

The specific angular momentum £ of a particle in a 
circular orbit with a radius r = xM in Schwarzschild 
coordinates around a mass M (where here and hence- 
forth we use geometrized units in which G = c = 1) 
is given by £ 2 = x 2 M 2 /(x — 3). The specific energy 
of this particle is E = (x — 2)/ y 'x{x — 3), but more 
generally for a particle of angular momentum £ in in- 
stantaneously azimuthal motion the specific energy is 
E=y/(1- 2/x)[l + £ 2 /{x 2 M 2 )\. 



Now suppose that a fraction e of the mass is taken 
away instantly. As before, we argue that the angular 
momentum is conserved, so that 

* 2 = x' 2 (l-e) 2 
x-3 x f -3 [ } 

where we define x' = r'/M, where r' is the radius after 
circularization and M is the final mass of the merged 
black hole. The specific energy immediately after mass 
loss can then be rewritten as 



E h 



l(x-2)[x-2(l-e)] 



x(x — 3) 



and the final specific energy after circularization is 

x'-2 



(5) 



(6) 



y/x'(x' ~ 3) 

These expressions include the rest-mass energy of the 
fluid, so to calculate the fractional change in binding en- 
ergy, we use 

(-Emit — -Efin) 



AE/E tot = 



(1 - E h 



(7) 



As an example of the degree of enhancement achieved 
relative to the Newtonian limit, consider a fractional 
mass loss e = 0.05. In the Newtonian limit, AE/E to t = 
0.0025. For x = 6 (the location of the ISCO), AE/E tot = 
0.0679, or 27 times larger. For x = 7, AE/E tot = 0.0295. 
The relative enhancement is even greater for smaller frac- 
tional mass losses. For e = 0.01, AE/E tot = 0.0001 in 
the Newtonian limit. In contrast, AE/E tot = 0.00582 
for x = 6 and AE/E tot = 0.00131 for x = 7. The en- 
ergy release is tremendously enhanced by proximity to 
the ISCO. 

The orbit-crossing timescales can also be reduced sim- 
ply because the gradient of the radial epicyclic frequency 
is much greater near the ISCO; indeed, at the ISCO itself, 
n = 0. We expect this speedup, as well as the enhance- 
ment of energy release, to apply for the real spacetime 
as well as for Schwarzschild, although of course the exact 
numbers will be different. 

Although these analytic expectations set the basic en- 
ergetic stage, numerical simulations are essential to eval- 
uate the details. In particular, the preceding analytic 
arguments did not address how fluid quantities, such 
as the pressure, would react to and modify the induced 
epicyclic motion and the eventual energy dissipation. In 
the next section we therefore present the results of our 
hydrodynamic and magnetohydrodynamic (MHD) simu- 
lations. We find that, as expected, in our models shocks 
and thus efficient thermalization of energy only occur for 
disks that are geometrically thin enough that the sound 
speed is less than the speed of radial epicyclic oscillation. 
This limits the accretion rate and also means that the 
disk viscous time is long. Therefore, the binary decouples 
from the disk when the disk is still far from the ISCO, 
so the relativistic energy enhancements do not apply. In 
addition, we find in our simulations that any luminosity 
enhancement is localized and minor. We do find, how- 
ever, that the increase in disk inner radius caused by the 
abrupt mass loss results in a large decrease in luminosity 
that is potentially detectable by future instruments. We 
focus on this effect in 8 4. 
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3. HYDRODYNAMIC AND MHD SIMULATIONS 

Our primary goal in conducting these simulations is 
to examine the extent to which shocks induced from a 
binary merger appreciably heat the circumbinary disk 
or generate other potentially observable electromagnetic 
signatures. It is reasonable to assume that any circumbi- 
nary disk will have evolved into a turbulent state, so we 
therefore begin by simu lating a baseline 3D MHD disk us- 
ing the methods (§ |3.1| ) and initial parameters (§ |3.2| ) de- 
scribed below. In the baseline model, accretion is driven 
naturally by turbulence a ssociated with the magn etoro- 
tational instability (MRI, Bal bus fc Hawley||1991 ), pro- 
ducing a disk that extends down to the inner regions of 
the potential near the IS CO. As discussed in the previous 
section, it is important to model this inner disk since that 
is where both the energy enhancement and timescale are 
optimal for generating a detectable merger signal. We 
then use the output of this baseline model to seed the 
initial conditions in our various mass -loss models, which 
are described in detail in Section [3T2I 



TABLE 1 

Simulations of Circumbinary Accretion Disks 
Experiencing a Central Mass Loss 



3.1. Numerical Methods 

To simulate the response of a circumbinary accretion 
disk to an instantaneous central mass loss event, we em- 
ploy the ZEUS-M P code (version 2), the de t ails of which 
are described in Stone & Norman (199 2a|b ), |Stone et al.| 
( [1992D , and |Hayes et al.| ( |2006p ~This code uses an Eu- 
lenan finite difference scheme to solve to second-order 
spatial accuracy the equations of ideal compressible fluid 
dynamics. We conduct our simulations using spherical 
coordinates (R, 6>, </>) in both full 3D MHD and "2.5D" 
hydrodynamic modes, the latter of which includes ax- 
isymmetric azimuthal velocities. The simulations feature 
a polytropic equation of state (p oc p 1 ) with 7 = 5/3. 
Additionally, we include a protection routine to impose 
a density floor p min at a value 10 -7 times that of the 
initial midplane density. We have also incorporated into 
th e code an Alfven spee d limiter of the type described 
in |Miller fc Stone] ( |2QQQ[ ). This algorithm mimics the ef- 
fect of the displacement current in the MHD simulations, 
naturally preventing the Alfven speed from exceeding the 
speed of light. 

Gravity in our simulations is handled through the in- 
clusion of a point mass at the computational grid origin, 
and self-gravity of the disk is not modeled. Since most 
of the gravitational wave emission in such a merger takes 
place in the last half-orbit, the binary is sufficiently close 
that it is reasonable to treat the pre-merger system as a 
single point mass. This approach does, however, preclude 
our simulations from experiencing any of the natural ef- 
fects of a true binary potential, such as the production 
of spiral density waves in the disk. We assume such phe- 
nomena will be of minor importance in the overall behav- 
ior of the system. We have modified the code to make 
the gravitati onal potential pseudo- Newt onian in form, as 
described by iPaczynski & Wiital (fl980|) : 



GM 



GM 



(8) 



This form of the potential features an IS CO at i? = 6M, 
the presence of which naturally augments the energy re- 
lease from a binary merger (see Section 2.3). Due to the 
relatively short timescale over which gravitational waves 



ID 


Type of Simulation 1 


Resolution 2 


e 3 


H10 


Hydrodynamic 


768 x 256 


0.10 


M10 


MHD 


768 x 256 x 32 


0.10 


H5 


Hydrodynamic 


768 x 256 


0.05 


HI 


Hydrodynamic 


768 x 256 


0.01 


Ml 


MHD 


768 x 256 x 32 


0.01 



1 All simulations are initialized from a single 3D MHD model. Mod- 
els M10 and Ml use the full data as initial conditions while models 
H10, H5, and HI input only the azimuthally averaged hydrodynamic 
variables. 2 Resolution is listed as R x 9 (x<p). The physical extent of 
the grid is R G (4M, 100M), 6 G (0.05tt, 0.95tt), G (0,tt/6). 3 We de- 
fine e as the fraction of the total initial mass lost through the binary 
merger, and we model the central mass loss as an instantaneous event. 
In the figures and discussion, the mass loss is described as taking place 
at a time t = 0. 

are emitted, the merger event itself is modeled as a dis- 
continuous change of the central point mass. 

All of our simulations span a radial domain R G 
(4M, 100M). This range is occupied by 768 zones for 
which AR increases logarithmically, with A.R m i n = 
0.02 M at the inner radial boundary. Both the inner 
and outer radial boundaries are used to enforce zero- 
gradient outflow. This is accomplished by setting fluid 
variables in the boundaries to their values in adjacent 
grid zones, with the additional condition that radial ve- 
locities are restricted so that flow on to the grid from the 
boundaries is not allowed. In the 6 direction, resolution 
is concentrated in the inner seven disk scale heights on 
either side of the midplane. There, AO w 2.1 x 10 -3 ra- 
dians, while AO increases logarithmically outside of this 
region so that the total grid spans G (0.057T, 0.957r) with 
256 zones. A reflecting boundary is employed in due 
to the proximity of the boundary to the coordinate pole, 
which would naturally reflect any approximately axisym- 
metric flow. This non-uniform grid structure has the ad- 
vantage of reducing the potentially unphysical influence 
of the outer grid boundaries by placing them physically 
far away from the regions of interest. For those simula- 
tions that are fully three-dimensional, <j) G (0,7r/6) with 
a uniform A<p « 0.016 radians. This should be sufficient 
to allow fully three-dimensional MHD behaviors while 
avoiding the excessive computational cost of spanning a 
full 2tt in azimuth. In both the axisymmetric and 3D 
simulations, periodic boundaries are enforced in the <p 
direction. 

3.2. Model Parameters and Initialization 

As mentioned above, we conduct a single 3D MHD 
baseline simulation to generate the initial conditions for 
all of our models. This simulation is run for ~ 6160 M, 
corresponding to ~ 100 Tisco? where Tisco is the orbital 
period at the ISCO. Since we are interested primarily 
in the disk properties just before and following binary 
merger, we define t = as the time of the mass-loss 
event. In this timescale, all data shown for reference at 
times t < reflect the output of the baseline simulation, 
while all data for which t > describe responses of our 
various disk models to the mass- loss event. 

The baseline simulation is initialized with complete az- 
imuthal symmetry and an aspect ratio h/r = 0.05, where 
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Fig. 1. — Snapshots of the midplane pressure taken at t = 384 M (solid line), where the central mass loss has taken place at t = (dashed 
line). Each of models H5, H10, and M10 shows a peak in pressure steep enough to represent a shock driven by the oscillation of the inner 
disk. In contrast, the low fraction of mass lost in model HI produces a much weaker and broader signal. 

r = RsinO is the cylindrical disk radius. The initial den- 
sity and pressure profiles are given by 

/ ^ ^ ( cos 2 \ 

p(R, 0) = po exp - - x9 ; , (9) 



2(/i/r) 2 



sin 



and 



p(R,6) 



GMR(h/r) 2 sin 2 
(^-2r g ) 2 



P(R,0), 



(10) 



where po is the initial density in the disk midplane 
(i.e., at = 7r/2). These forms result in an initial disk 
of constant midplane density and a radially decreasing 
temperature profile. The initial disk is also in approxi- 
mate hydrostatic equilibrium, although this is disrupted 
by the development of turbulence driven by the MRI. 
The initial disk velocity is entirely azimuthal with 



V<f> 



VgWt 

r — 2rv 



vr = v = 0. 



(11) 



The initial magnetic field configuration consists of a 
series of field loops, each of which is contained within 
an R — plane. This field was derived from a vector 
potential of the form 



A p 



1/2 



sin 



2nr 

dloopO) 



[dioop(r)] 2 (12) 



where Aq is a scaling constant chosen such that the ini- 
tial maximum magnetic field strength corresponds to a 
ratio of gas-to-magnetic pressure of f3 ~ 3000. The sim- 
ple function d\ oop = 1.5 + r/20 is included to increase the 
width of the field loops with increasing r. The develop- 
ment of the MRI should be reasonably insensitive to the 
exact details of this field except insofar as the field has no 
net flux over scales larger than a few gravitational radii 
in the inner disk. By the time this baseline simulation 
has run for 100 Tisco? the inner third of the accretion 
disk (i.e., out to ~ 35 M) has evolved to a turbulent 
state through the action of the MRI. It is this state that 
is imported as initial conditions into our various mod- 
els of mass loss from a binary merger. For those fully 
3D MHD models (denoted with an "M"), we use exactly 
the output of the baseline simulation as our initial condi- 
tions. The "2.5D" hydrodynamic (denoted with an "H") 
models, on the other hand, incorporate the azimuthally 
averaged output from the baseline simulation, neglect- 
ing the magnetic field information. While this initial 
state is slightly out of equilibrium because of the omis- 
sion of magnetic and/or non-axisymmetric stresses, the 
gross similarity between the behaviors of the MHD and 
hydrodynamic models suggests that this is not a signifi- 
cant problem. The only other parameter varied between 
models is the amount of central mass lost to gravitational 
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Fig. 2. — Snapshots of the density in model H10 at various times. The colorbar is logarithmic, ranging from 10 -4 po (dark) to po (light). 
The mass-loss event takes place at t = 0. The initial wave from the mass-loss event is visible at t = 134 M while the stronger shock is 
visible at t = 384 M and t = 1384 M. Features such as the shocks and rarefactions propagate radially outward in the disk (from left to 
right in these images). 



waves, referred to in the model names as a percentage 
(e.g., model "M10" is an MHD disk that loses 10% of 
its central mass). Each of these simulations is run for at 
least 20 IS CO orbital periods after the mass- loss event 
has taken place. The simulations are summarized briefly 
in Table 1. 

3.3. Simulation Results 

We now describe the results of our simulations of accre- 
tion disks surrounding central objects that have lost mass 
through the prompt emission of gravitational waves. It is 
worth noting that we will frequently use descriptors such 
as the gravitational radius and the ISCO orbital period, 
both of which are defined in terms of the central mass. 



Rather than having to adjust our spatial and temporal 
reference scales after the mass- loss event, all distances 
and times are given in units of the initial central mass. 

3.3.1. Disk Response and the Formation of Shocks 

As described earlier, the basic response of any orbiting 
fluid element to a sudden reduction in the central mass 
is radial oscillation. Since the disk material is suddenly 
at the pericenter of what will become its new orbit, the 
initial radial motion after the mass-loss event is directed 
outward (v r > 0). As expected, all of our simulated disks 
feature this generic behavior, and the amount of outward 
disk motion is observed to depend upon the fraction of 
central mass that has been lost. Since the hydrodynamic 
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Fig. 3. — Comparison of the mass accretion rate across the ISCO for hydrodynamic (left) and MHD (right) models. Mass-loss values of 
10% (solid lines), 5% (dashed line), and 1% (dot-dashed lines) are shown, with the mass-loss event taking place at t = 0. In hydrodynamic 
disks, accretion never resumes after the mass-loss event. In model M10, however, there is a initial drop in accretion rate, followed by a 
gradual resumption of accretion over the viscous timescale. The initial accretion rate prior to the binary merger is shown at times < for 
reference. 



"H" models and MHD "M" models feature very similar 
large-scale disk dynamics, we will focus in this section 
primarily on the simpler hydrodynamic models, noting 
differences in the MHD cases at the end of the section. 

We first discuss model H10, which features the largest 
amount of central mass loss and therefore the most dra- 
matic disk response. Defining disk material by a density 
cutoff of Pdisk > O.Olpo? the inner edge of the H10 disk 
is seen to migrate out from r ^ 4 M to r ~ 12 M in 
just over one ISCO orbit. The progress of the inner edge 
of material then stalls out at this radius while the ini- 
tial perturbation from its motion continues to propagate 
radially outward. The innermost edge then returns radi- 
ally inward, reaching a minimum radius of approximately 
r ~ 8 M before beginning the next cycle of oscillation. 
This first oscillation is completed in a total of roughly 
four ISCO orbits before the next cycle begins. We should 
mention that the entire body of the disk is undergoing 
similar oscillations, but that the inner disk oscillation 
timescale is shortest and thus most easily identified. 

Each of the first two inner disk oscillations in model 
H10 drives a signal radially outward through the body 
of the disk. The detailed structures of these features 
change as they propagate through the disk, but they are 
generally identifiable at later times. The lower- left panel 
of Figure [I] shows the instantaneous midplane pressure 
profile for model H10 at t = 384 M, approximately six 
ISCO orbits after the mass- loss event. The initial disk 
perturbation is now seen as a small peak near r = 28 M, 
while the second signal forms a very prominent peak at 
r = 12 M. The two signals are separated by a rarefaction 
with a pressure nearly an order of magnitude below the 
pre-loss pressure value (shown in Figure fl] as a dashed 
line) . This rarefaction, caused by the local depletion of 
disk material into the regions of enhanced density and 
pressure, also propagates radially outward, bracketed ra- 
dially by the two pressure peaks. The steepness of the 
leading edge of the inner pressure feature in this figure 
suggests that this oscillation is driving a shock radially 
outward into the disk. Similarly, Figure [2] shows a ren- 
dering of the disk density at four different times during 
the evolution of model H10. In this figure, the initial per- 



turbation is only easily seen at t = 134 M, but the second 
stronger perturbation is identifiable at both t = 384 M 
and t = 1384 M. We also note that the disk material 
behind the second density /pressure enhancement can be 
seen expanding vertically out of the midplane, suggest- 
ing that heating is taking place. As seen in the bottom 
two panels of Figure [2| this heated material actually gets 
ahead of the midplane shock front by pushing quickly 
through the relatively diffuse material at higher latitudes 
before spiraling back into the disk. After the first two 
oscillations, subsequent individual oscillation events be- 
come much more difficult to identify (see the bottom- 
right panel of Figure |2| , and the inner edge of the disk 
gradually settles down to a radius r ~ 9 M. 

In contrast, Model HI does not feature enough central 
mass loss to form shocks. The upper-left panel of Figure 
[I] shows the pressure profile of model HI, also approxi- 
mately six ISCO orbits after the mass loss. There is an 
identifiable feature near r ~ 15 M that corresponds to the 
second disk oscillation, but it is too broad to be a shock. 
Moreover, it is difficult to identify the initial pressure 
signal, which no longer stands out above the background 
profile. The inner edge of the HI disk moves at most a 
few gravitational radii as a result of the mass- loss event, 
and this motion is simply insufficient to drive shocks into 
the disk. Unsurprisingly, Model H5 features behavior 
that is bracketed by the other two models. Specifically, 
the upper-right panel of Figure [T] shows that H5 appears 
to form a shock on the second disk oscillation, but also 
that the initial pressure signal does not stand out. The 
rarefaction between these two peaks is also visible, but 
reduced from that of H10. 

Combined, these three simple models suggest that 
there is a lower limit to the amount of central mass loss 
capable of generating shocks in an accretion disk. This 
limit is not surprising, given that the requirement for 
a shock to form is v r > c s , assuming that v r is also an 
approximation of relative velocity in the system. Rewrit- 
ing c s ~ (h/r)V(p and noting that v x jv^ w e for mass-loss 
events, the requirement becomes e > h/r. Although we 
have not fully explored the parameter space of mass-loss 
values, our simulations are consistent with this limit in 
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Fig. 4. — Comparison of the bremsstrahlung luminosity as measured between R = 5 — 35 M for hydrodynamic (left) and MHD (right) 
models. Mass-loss values of 10% (solid lines), 5% (dashed line), and 1% (dot-dashed lines) are shown, with the mass-loss event taking place 
at t = 0. In both hydrodynamic and MHD disks, there is barely any increase in the luminosity. There is, however, a drop in luminosity 
that depends upon the fractional mass lost, followed by a gradual recovery. The initial luminosity prior to the binary merger is shown at 
times < for reference. 



that shocks only fail to form in the one model for which 
e < h/r, and we expect this simple trend to apply more 
generally to any thin disk aspect ratio. As we will dis- 
cuss in Section 4.1, this limits us to considering shocks 
only in low- luminosity systems for which h/r < e, which 
is itself much less than unity. 

Finally, we point out that the 3D MHD simulations 
feature much of the same behavior as the hydrodynamic 
models. Figure [T] illustrates that models M10 and H10, 
for example, produce very similar pressure profiles at 
comparable times. Specifically, model M10 clearly fea- 
tures a shock and an identifiable initial peak very simi- 
lar to those of H10. Likewise, the behavior of the inner 
disk edge is quite similar for the first few oscillations. 
Model Ml is similarly akin to model HI in that neither 
forms shocks. In fact, the similarity of hydrodynamic and 
MHD simulations at early times is good evidence that the 
azimuthally averaged, non- magnetized initial conditions 
used in the hydrodynamic models did not induce stresses 
that were significant compared to those triggered by the 
mass loss. There is, however, one notable difference be- 
tween hydrodynamic and MHD models that reveals itself 
at later times. Whereas the inner edge of model H10 set- 
tles down to r ~ 9 M, for example, the inner edge of 
M10 gradually begins to migrate inward. This is due 
to the resumption of accretion, which is accomplished in 
M10 through turbulence driven by the MRI. Although 
model H10 features an initially turbulent disk profile, the 
usual means of driving turbulence are not included in the 
physics of the hydrodynamic models. While this differ- 
ence does not strongly affect the formation of shocks, it 
does affect possible observational signatures of mass loss, 
as we will discuss in the following section. 

3.3.2. Simulated Observations of Mass Loss 

While our simulations do not include the physics 
needed to assemble rigorous models of radiation produc- 
tion and transfer in accretion disks, we can construct 
simple proxies for real light curves. Here, we consider 
the mass accretion rates and bremsstrahlung luminosi- 
ties derived from our model disks to explore how these 
quantities reflect the central mass loss and disk dynamics 



described in the previous section. 
First, we construct the mass accretion rate as 



M 



-L 



R(-v R )pdS, 



(13) 



ISCO 



where the integral is taken over the surface formed by 
the radius of the ISCO. Figure [3] shows this quantity 
plotted for all five disk models. As described in the pre- 
vious section, the hydrodynamic models lack the physics 
needed to properly resume accretion after the mass-loss 
event. The oscillation of the inner edge of model HI 
intersects the ISCO, creating a pattern of semi-regular 
fluctuations, but only a negligible fraction of this inflow 
is accreted onto the compact object. The MHD models, 
in contrast, show much more interesting and realistic ac- 
cretion behaviors. Both models Ml and M10 feature a 
sudden depression in M just after the mass- loss event. 
This represents a decrease of over two orders of magni- 
tude in M as accretion is temporarily disrupted. In the 
case of model Ml, this reduced accretion rate persists 
for at most two ISCO orbits (~ 120 M) before accretion 
across the ISCO rapidly resumes. Model M10, however, 
shows a much longer-lived decrement in mass accretion 
rate. After the initial drop, the accretion rate in M10 
grows gradually at a rate that, if extrapolated, would 
match the initial accretion rate after a time ~ 2000 M. 
This is in fact comparable to the viscous timescale of 
the inner disk, as measured after the mass loss has taken 
place. 

We also examine the bremsstrahlung luminosity, which 
we compute from 



L 



brem 



(14) 



where ebrem = p 2 T x l 2 is a measure of the total 
bremsstrahlung emissivity. This is admittedly a some- 
what naive treatment of the complex radiative processes 
expected in real disks, particularly those in the opti- 
cally thick regime. We argue, however, that the effect 
of opacity would only be to dilute or delay a radia- 
tive signature of mass loss, and that the optically thin 
regime constitutes an important test case. Specifically, 
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Fig. 5. — Comparison of the bremsstrahlung luminosity as measured between R = 20 — 25 M for hydrodynamic (left) and MHD (right) 
models. Mass-loss values of 10% (solid lines), 5% (dashed line), and 1% (dot-dashed lines) are shown, with the mass-loss event taking 
place at t = 0. In both hydrodynamic and MHD disks, there is an initial increase in luminosity, followed by a drop that depends upon the 
fractional mass lost. This is then followed by a sharp recovery as the shock passes through, followed by another drop. Although there are 
short-lived increases in the luminosity of this single annulus, Figure [4] shows that such features are not visible when the luminosities are 
summed over a larger region. The initial luminosity prior to the binary merger is shown at times < for reference. 



if our optically thin models fail to produce a signifi- 
cant bremsstrahlung signal, further radiative processing 
is irrelevant. If, on the other hand, the bremsstrahlung 
emission profile proves to be a clear indicator of mass 
lost through a binary merger, then we must include the 
caveat that such emission will be subject to radiative re- 
processing. We will discuss this issue further in Section 
4, where we examine in detail the observability of signals 
from mass loss. 

Figure [4] shows for all models the bremsstrahlung lu- 
minosity as constructed for a thick spherical shell rang- 
ing from R = 5 — 35 M. The most interesting aspect 
of this figure is that none of the models feature a pro- 
nounced increase in total bremsstrahlung emission in this 
region. Naively, one might have assumed that the ac- 
tion of shocks passing through this region would only 
increase the luminosity. Instead, however, we see that 
the reduction in pressure and density in the rarefaction 
(visible in Figures fl] and J2| actually works to reduce 
the total bremsstranlung luminosity, followed by a re- 
covery as the shock passes through the region. Simulta- 
neously, the oscillations deplete material from the inner 
regions of our integration volume, also acting to lower 
the luminosity. The net effect of these processes is a to- 
tal bremsstrahlung luminosity that essentially never in- 
creases above the baseline level. 

In order to examine only the effects of shocks and 
rarefactions on this disk, we show in Figure [5] the 
bremsstrahlung emission in a narrow annulus ranging 
from R = 20 — 25 M, well outside of the range of mo- 
tion of the inner edge of the disk. In fact, Figure [5] clearly 
shows indications of each of the pressure and density fea- 
tures seen in Figures [I] and |2j Considering for the mo- 
ment models M10/H10, the luminosity increases at first 
by about 10 — 20% as the perturbation moves through the 
volume. This is followed by a sharp J> 90% net drop in 
luminosity as the rarefaction passes through. The next 
compression causes the next luminosity increase, which 
can be significant for models H10, H5, and M10, all of 
which produce shocks. We see, however, that shocks 
only increase the luminosity by at most ~ 40% above 



its original value. Furthermore, this enhancement is not 
long-lived. Once the next rarefaction arrives, the lumi- 
nosity drops again, and it is not obvious that the net lu- 
minosity will be significantly increased over long times. 
While Figure [5] confirms that the shocks and rarefac- 
tions are behaving as expected, we emphasize that the 
observed trends in luminosity are essentially local phe- 
nomena. When the emission from several such annuli is 
summed, as shown in Figure [4] for example, we see no 
significant increase in bremsstrahlung luminosity. This 
suggests that any local increases in luminosity are con- 
stantly being offset by decreases elsewhere in the disk. 

In fact, the drop in bremsstrahlung luminosity is the 
dominant feature in Figures [4] and [5] As mentioned, the 
M10/H10 models in particular feature a J> 90% reduc- 
tion in luminosity, followed by a gradual recovery. This 
effect is much more dramatic than the minor luminos- 
ity increase produced by the initial perturbation, as seen 
in Figure [5] Moreover, we see that the presence of this 
luminosity dip does not require shock formation. Specif- 
ically, the Hi/Ml models feature a ~ 20% reduction in 
luminosity and do not form shocks. This is significant 
since models without shocks are no longer constrained 
by the e > h/r requirement and may thus be related to 
more luminous systems. 

A comparison of Figures [4] and [5] also illustrates that 
the recovery time of the prompt bremsstrahlung luminos- 
ity drop is determined at least in part by the size of the 
region we choose to examine. The bremsstrahlung lumi- 
nosity profile of M10 in Figure [4j for example, resembles 
the profile of M, but the reasons for these similar profiles 
are quite different. In the case of M, the recovery time 
was determined by the physical parameters that feed into 
the viscous timescale. In the case of the bremsstrahlung 
luminosity, Figures [4] and [5] show that the recovery time 
depends upon the size of the sampling volume. From an 
observational point of view, this implies that, in the case 
in which the observed luminosity is predominantly pow- 
ered by accretion onto the binary, one may be able to 
measure the viscous timescale of the circumbinary disk 
from the light curve and thus constrain the structure of 



10 



the disk. Alternately, if the true luminosity is dominated 
by emission from the disk itself (here estimated in the 
form of bremsstrahlung emission), the drop in flux may 
be detected as a consequence of the mass loss. The length 
of the drop in that case, however, would not reflect the 
viscous timescale. 

4. OBSERVABILITY 

The observability of signals from mass loss due to grav- 
itational wave emission depends on the flux and duration 
of the event, but also on the angular localization possible 
from LISA observations. In this section we explore var- 
ious possible signals and observing strategies. We focus 
primarily on the energetics (as opposed to detailed spec- 
tral modeling), but note that the strongest signals are 
likely to be found in the X-ray band or, in some scenar- 
ios, the radio band. Both of these bands are dominated 
by processes associated with the innermost parts of the 
accretion disk, where the effects of the abrupt mass loss 
are the strongest. Furthermore, they can penetrate sig- 
nificant columns of circumnuclear matter that may be 
present. On the other hand, ultraviolet, optical and in- 
frared variations are likely to be relatively weak due to 
dilution by the outer accretion disk and starlight and will 
be much more susceptible to absorption. Note however 
that the precursor signature in any band may be difficult 
to observe in Compton-thick AGN due to obscuration. 

We briefly comment upon the mechanism by which dis- 
turbances in the accretion disk may be reflected in the 
observed X-ray and radio fluxes. X-rays are likely pro- 
duced via thermal Comptonization of softer disk photons 
by a corona of energetic electrons. The energy source 
for this corona is believed to be the accretion disk it- 
self; abrupt changes in the energetics of the disk are very 
likely to lead to abrupt changes in coronal heating and 
X-ray emission. Radio emission is associated with jets 
formed in the inner accretion disk or, possibly, the black 
hole magnetosphere. Even in the black hole magneto- 
sphere scenario where the p ower source of the jet is ul - 



timately the black hole spin ( |Blandford fc Znajek|[l977 ), 
the magnetic fields are generated by currents in tne disk 
and hence the jet (and associated high-frequency radio 
emission) will quickly respond to changes in the inner 
accretion disk. Thus, both X-ray and radio emission are 
expected to drop significantly if the inner accretion disk 
is disrupted or recedes and should increase as the inner 
accretion disk is replenished. 

As a practical matter, it is very important to recog- 
nize that AGN are intrinsically variable and have many 
different spectral distributions. For example, the narrow- 
line Seyfert 1 galaxy NGC 4051 exhibits rms variation of 
tens of percent in X-ray emission from ~ 10 -8 — 10 ~ 3 Hz 
( |McHardy eTH] [2004] ) . AGN light curves are also 
known to b e non-Gaussian, sometim es showing large 
flares (e.g., [Leighly fc 0'Brien|[l997l ). Therefore, mere 
detect ability of the flux excess or deficit from a given ef- 
fect is not sufficient to make the source stand out from 
other AGNs. If, however, a predicted flux change hap- 
pens at the point of merger as determined by gravita- 
tional wave observations, the probability that this is an 
unrelated source drops dramatically. In addition, such 
simultaneity between electromagnetic and gravitational 
radiation would constrain strongly any deviation of the 
speed of gravitational waves from the speed of light. 



We now explore the observability of extra emis- 
sion from circularization and dissipation, and from the 
prompt drop in emission caused by mass loss. 

4.1. Observation of energy from circularization 

The extra energy release due to circularization of ellip- 
tical orbits is not likely to be observable. The first reason 
is that our analytic treatment shows that the additional 
luminosity due to the circularization process is at most 
10% of the accretion luminosity. Indeed, our MHD simu- 
lations show that the maximum bremsstrahlung emissiv- 
ity only comes back to the level achieved prior to mass 
loss, and that there is not a series of peaks that could be 
picked out easily from the normal MHD turbulence. Ad- 
ditionally, this is independent of any radiative reprocess- 
ing that would presumably only further muddle signals 
of interest. 

The second reason is that for shocks to develop and 
thus for energy to be thermalized efficiently, our simula- 
tions show that the fractional disk thickness needs to be 
h/r < e. This limits the radius of the inner edge of the 
disk. To compute this we follow Armitage fc N atarajan 
( |2002 ), who note that the viscous rate at which the inner 
edge of a Shakura-Sunyaev disk moves inwards is 



avK 



(15) 



where a is the viscosity coefficient and vk = ^jGMja 
at a semimajor axis a, whereas the rate at which the 
semimajor axis of a circular binary shrinks is 
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If we define a dimensionless parameter x by a = xGM/c 2 
and set h/r = 0.05(4?7) (i.e., suppose the disk is 
marginally capable of forming circularization shocks) 
then the viscous and binary shrinkage speeds become 



^visc = -0.02477 2 (a/0.4)x" 
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Setting the two equal gives 

x min « 20(4r ? )-°- 4 (a/0.4)-°- 4 



(17) 



(18) 



In these equations we scale by a viscosity parameter 
a = 0.4, which is the maximu m of the ran g e infer red from 
observations summarized by 
fore, in units of gravitational 



King et al.N2007 ). There- 
radii, the closest the disk 



could get to the merging binary is about 20 M, and this 
occurs for equal-mass black holes. In reality, the in- 
ner edge will be farther away, because for comparable- 
mass black holes the disk is truncated at about dou- 
ble the binary semimajor axis. This implies a mini- 
mum of about 40 M, or a bit further because at that 
larger radius, d v i SC is reduced from the value we calcu- 
lated above. Offsetting this somewhat is that the disk 
will drift in while the binary spirals in, so we will take 
dinner, min = 30 M as an estimate of how close the in- 
ner edge of the disk can be at merger and still produce 
shocks after mass loss. This implies a radiation efficiency 
of E = 1 - \u t \ = 1 - (x - 2)/y/x(x-3) = 0.016c 2 for 
a Schwarzschild spacetime, which should be a good ap- 
proximation at this distance. 
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To estimate the maximum mass accret i on rat e we use remains a poor candidate for detectability. 



the disk solution of Shakura & Sunyaev (1973). In the 
inner radiation-dominated portion of their solution (valid 
for all the radii of interest here), they find a disk half- 
thickness of 

h = 3.2 x 10 6 cm (M/l M Q )(M/M E )[1 - (r/6M)~ 1/2 ] 

(19) 

where M E = L E /(0.06c 2 ) and L E = 1.3 x 
10 38 (M/Mq) erg s _1 is the Eddington luminosity. In 
this region the half-thickness depends only weakly on ra- 
dius. At r = 30M, therefore, h/r w 0A(M/M E ). For 
e ~ 0.05 the requirement h/r < e means M/M E < 0.125. 
The maximum disk luminosity is therefore 0.016c 2 M = 
0.125 x (0.016/0.06)^ = 0.03L E . 

To be specific, consider an equal-mass binary of total 
mass 2 x 10 6 M accreting at this rate, at a redshift of 
z = 1 and thus a luminosity distance of d E = 6.6 Gpc (as- 
suming a flat universe with matter parameter Q m = 0.27 
and Hubble parameter Ho = .71). The maximum ob- 
served flux is then F max = 1.5 x 10 -15 erg cm -2 s _1 . 
According to our earlier estimate of Newtonian shocks, 
the flux variation due to mass loss would be a factor of 
~ 10 lower than this. The flux variation would there- 
fore be just a fewxlO -16 erg cm -2 s _1 , so even if all 
of the radiation emerges in X-rays it is only compara- 
ble to the dimmer sources i n the lMsec Chandra Deep 
Field (Giacconi et al. 2002). The duration of this ef- 
fect would only be roughly the orbital time at 30 M, or 
2tt x 30 3 / 2 « 1000 M, which for M = 2 x 10 6 M is only 
about 10 4 seconds. 

The most capable X-ray observatory currently planned 
that may operate contemporaneously with LISA is the 
International X-ray Observatory (IXO). However, even 
for IXO (and temporarily putting aside issues of angular 
localization and field-of-view; see §4.3), the flux limit of 
a 10 4 s integration is ~ 10 -15 erg cm -2 s _1 . So even 
if this one-orbit variation were distinguishable from the 
natural fluctuations caused by MHD turbulence, it is too 
faint to be practically detectable. 

Although we have focused on the case in which h/r ~ e, 
we should mention that smaller h/r ratios will generally 
allow for the development of more powerful shocks. Un- 
fortunately, such disks also have correspondingly larger 
inner edge radii and lower characteristic luminosities, so 
the prospects for detection are not necessarily improved. 
Finally, there is the strong possibility that accretion onto 
the binary system will not be perfectl y steady . As dis- 
cussed by, for exam pl e, I vanov et al. (|1999|), Arm itage 



cussed by, tor example, ivanov et al. (ryyy), Armitage 
fc Natarajanl (|2QQ2|), |m osavljev iT"fc R[mmey[^2 005 ) 



MacFadyen & Milosavljevic (2008|, and |Haiman et al. 



( 2009]), torques induced by the binary are likely to mod- 
ify the structure of the inner disk. While the dominant 
effect is expected to be a clearing out of disk material 
within approximately twice the binary orbital radius, as 
we have assumed above, it is also likely that such torques 
will modify th e density structure of the inner dis k and, 
as discussed in [Armita ge & N atarajan] ( 2002 ) and |Mac-| 
Fadyen & Milosavljevic (2008), drive spiral waves into 



the disk body. Although we have not explicitly modeled 
such features, it seems unlikely that they would achieve a 
signal of the energy release from circularization stronger 
than that which occurs in a steady disk. This particu- 
lar electromagnetic signature of binary merger therefore 



4.2. Observation of prompt drop influx 

In contrast, our simulations show that the drop in mass 
accretion rate caused by the outward movement of the 
inner edge of the disk can last for the viscous time at 
that radius. Because this effect does not require forma- 
tion of shocks, the h/r < e requirement can be relaxed. 
Therefore, the abrupt reduction in mass accretion could 
be observable in the form of a prompt drop in flux even 
in systems close to the Eddington luminosity. In these 
systems the viscous time is short enough that the inner 
edge of the disk can keep up with the binary all the way 
to the point of dynamical instability (note that if we set 
h/r ~ 1 and a = 0.1 and solve as in equation (17), we 
find that the disk stays with the binary until r 5 6 M). 
This increases the maximum flux for systems by a factor 
of 30, or up to F max ~ 5 x 10 -14 erg cm -2 s _1 , com- 
pared to the flux estimate for a geometrically thin disk 
in § 4.1. The recession of the disk edge following the 
black hole merger (and subsequent mass loss) will lead 
to order unity decrease in the total flux, i.e., a transient 
that is ~ 200 times stronger than the luminosity effect 
expected from circularization. At such high accretion 
rates, the viscous time is t v i SC ~ t rb/<^> implying that 
the flux deficit will last for 600 — 2000 M depending on 
the location of the inner edge. For a total binary mass 
of ~ 2 x 10 6 M and a redshift of z = 1 this translates 
to ~ 1 — 4 x 10 4 seconds. 

Of course, since we are now considering the cessation 
and then resumption of accretion, we need to include the 
bolometric correction factor linking this total flux with 
that in the X-ray band. For these small mass systems, 
we will assume that 10% of the total flux emerges in the 
X-ray band. Thus, the signal we seek is the following: an 
AGN with an X-ray flux of F x ~ 5 x 10 -15 erg cm -2 s _1 
which abruptly turns off at the moment of merger and 
then re-emerges on timescales of ~ 1 — 4 x 10 4 seconds. 

Localization of the source on the sky prior to the co- 
alescence of a supermassive black hole binary can be 
achieved if there is a clear electromagnetic signature as- 
sociated with the host galaxy. Examples of precursor 
signatures previously discussed in the literature include 
the disturbed morphology of the host galaxy and star- 
burst, as well as the precursor accretion episodes, and 
are typically expected to occur in case of binaries heav- 
ier than ~ 1O 7 M ( |Armitage fc Natarajanl|20Q2| |Dotti| 
et al.|2006 ). However, less-massive binaries (with masses 
< 5 x 1O 6 M0) are not expected to produce strong pre- 
cursor signatures, for several reasons. First, in the time 
required for binary black holes to coalesce (which could 
be up to a few Gyr), any gas trapped within the binary 
orbit would have been consumed, reducing the probabil- 
ity that the mass accretion rate is close to Eddington. 
Second, in this time it is also expected that the star- 
burst produced by the galactic merger would have run 
its course, and in addition it is unclear whether there 
are any remaining observable signs of a disturbed galac- 
tic morphology from the merger. Low-mass binaries are 
therefore expected instead to be followed by an electro- 
magnetic afterglow rising on the viscous time scale of the 



circumbina ry disk ( jMilosavljevic fc Phinney||2005 Dotti 
eraI1[2006l ) 



We do note that even low-mass binary systems may be 
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identifiable prior to coalescence and mass loss due to the 
increasing radiative efficiency as the inner radius of the 
disk decreases while the binary separation decays. For 
example, consider an equal-mass binary at a separation 
of 10M, where M is the total mass of the binary. Nu- 
merical simulations and simple estimates show that the 
inspi ral time from thi s point is approximately 1000M 
(e.g., Ba ker et al.|2QQ7 ). At lowest order th e inspiral tim e 
scales as the fourth power of the separation ( |Peters 1964[ ) , 
hence a binary with a separation of 20M would spiral in 
within 1.6 x 10 4 M, which for a M = 2 x 10 6 M binary 
at redshift z = 1 would appear to us as about four days. 
If the inner edge remains at twice the binary semimajor 
axis for the duration of the inspiral, this implies that the 
radiative efficiency and hence the luminosity may roughly 
double in that time. This is potentially observable, al- 
though we again caution that AGN are known to have 
flux variat ions of several tens o f percent over hours to 
days (e.g., McHardy et al.||2004 ). We are also assuming 
that the luminosity from accretion onto the individual 
black holes is small compared to the luminosity from the 
circumbinary disk. This is expected in the context of 
our model, where the density of the gas in the innermost 
region of the disk is very low in comparison to the disk 
mid-plane and the disk is assumed to be axisymmetric. 

In a realistic scenario the circumbinary disk will not 
be perfectly axisymmetric, and the black holes will not 
have equal mass. It is therefore possible that the rate at 
which the gas flows from the inner rim of the disk into 
the low density region around the binary may be signif- 
icant. Some of this gas w ill be accreted by the black 



holes ( |Hayasaki et al.||2008| ), while the remainder devel- 
ops eccentric orbits, and may collide with the inner rim of 
the circumbinary disk or l eave the binary system in form 



of high velocity outflows (lA xmitage & Natarajan 2002 



MacFadyen & Milosavljevic 2008). Accretion on one or 
both members of the binary and gas stream collisions 
with the disk may give rise to quasi-periodic outbursts 
on the orbital time scale. In the case of the binary ex- 
ample considered above this translates into a variability 
time scale of about ~ 3 hr in the frame of the observer. 
The variability time scale is expected to decrease mono- 
tonically as the orbit of the binary decays. In order for 
quasi-periodic outbursts to stand out above the circumbi- 
nary disk emission, their luminosity would have to be at 
least at the level of F m ax- Moreover, if a correlation can 
be demonstrated between pre-coalescence quasi-periodic 
outbursts and gravitational wave emission, this would be 
a telling sign that th e electromagnetic s ignature is related 
to the binary orbit ( |Kocsis et al.||2008 ). 

We emphasize that the observation of a drop in mass 
accretion rate, unlike the situation described in Section 
4.1, assumes the maintenance of accretion during the in- 
spiral period. If the inner edge of the accretion disk is 
stalled at large radii (r ~ 40 — 1 00M), as is expected for 
geometrically thin disks (e.g., [Armi tage fc Natarajan 



2002; |Milosavljevic fc Phirme y 2005), then the associ- 



ated dropoff in accretion would not directly reflect the 
merger event. As we have noted, however, systems that 
accrete near the Eddington rate should naturally feature 
disks with aspect ratios of order unity. Such disks should 
follow the binary system to small radii and, in the ab- 
sence of appreciable torques, undergo accretion on to the 
central object (s). If the binary torques on the disk are 



sufficient to substantially retard accretion, however, the 
associated accretion behavior would again fail to reflect 
the merger itself. It is only for the case in which accre- 
tion is modified but not completely prevented, such as 
the quasi-periodic outburst case described above, that a 
torqued disk could feature accretion behaviors directly 
indicative of a binary merger event. 

On the other hand, a dropoff in emission from the body 
of the disk, such as that discussed in Section 3.3.2, re- 
quires no constraining assumptions about accretion. In 
that section we used bremsstrahlung as a simple, generic 
proxy for disk emission and found that the luminosity of 
the inner r ~ 30M of the disk was reduced by a factor of 
2 (5) for e = 0.05 (0.10). In this scenario, an important 
caveat is that the inner portions of real disks are probably 
optically thick and would therefore reprocess this radia- 
tion. To address this issue, we estimate the approximate 
radiative diffusion time tdiff ~ hr/c, where h is again the 
disk half-thickness and r is the optical depth, in this case 
due to Thomson scat tering. Again applying the model of 
|Shakura fc Sunyaev| ( |1973| ) for the radiation-dominated 
inner portion of an accretion disk, 



1.8 (r/6M) 3 / 2 



a(M/M E )[l - (r/6M) -1 / 2 



(20) 



Combining this with the expression for disk height given 
in (19), we find that 



t dm ~ 4.4 x 10 4 s ( — 

a 



M 



2 x 10 6 M, 



r \3/2 



Gom) 



(21) 

For our fiducial values of a = 0.1, M = 2 x lO b M , and 
r = 30M, tdiff ~ 12 hours, which, at a redshift of z = 1, 
would appear to us to take a day. This is only a few times 
the local orbital period and rapid enough that scattering 
should not significantly inhibit the prompt emission of 
radiation from the disk body. 

We note that the example flux estimates in this section 
were calculated for a binary of mass 2 x 10 6 M . More 
massive binaries, however, would coincide with more lu- 
minous EM sources and may thus be easier to localize on 
the sky. In addition, the flux deficit in the light curves 
of these objects would last for a proportionally longer 
time. If such systems are too much more massive than 
our fiducial value, however, they will begin to fall outside 
of the LISA detection band. 

Although the uncertainties in the sub-parsec structure 
of a circumbinary region are reflected in the uncertain- 
ties of the associated signatures, it is plausible that the 
EM signatures are bracketed by the two simple scenarios 
described above. We therefore expect that the prompt 
drop in flux due to the mass loss is detectable only if 
the proper binary system can be identified. To deter- 
mine the likelihood of this, we now consider the proba- 
ble error boxes for such systems and requirements for the 
electromagnetic sky surveys of the future. 

4.3. Angular localization of coalescences 

A number of researchers have estimated the angular 
uncert ainty of black hole coalescences from LISA data 
{e.g., |Cutler| [19981 ISintes fc Vecchiol |2000l |Moore fc 



Hellings 2002; Hughes 2002; Vecchio 2004; Berti et al 



2005| |Holz fc Hughes 2005| |Lang fc Hu ghes 2006, 2008). 
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The size of the final error region, after inspiral, merger, 
and ringdown, will depend on the direction, redshift, 
black hole masses, and other factors, but also on whether 
the spin axes of the black holes are aligned with each 
other and with the o rbital axis. If they are not, then 



Lang & Hughes ( 2QQ6| ) show that precession in the final 
day of coalescence between typical supermassive black 
holes can reduce the major and minor axes of the error 
ellipses by factors of a few each. With this reduction the 
median major axis at redshift z = 1 for primary masses 
m the - 10 5 -10 7 M range is ~ 15'— 45' and the median 
minor axis is ~ 5' — 20' . However, models of supermas- 
sive black holes in gas-rich environments (i.e., the ones 
most favorable for high accretion rates and detectability) 
suggest that when the holes are still hundreds of parsecs 
from each other and thus their dynamics are dominated 
by the gas instead of each other, gas torques will align 
the spins with each other and with the orbital axis (Bog- 



danovic et al. 2007). If the axes are not aligned then 
the final error box near the point of merger is typically 
a few hundredths of a square degree at redshift z = 1. If 
the axes are aligned, then the lack of pre cession degrades 



this t o a few tenths of a square degree (Lang & Hughes 
2006]). 

Whether or not the axes are aligned, a day before 
merger the error region is likely to be in the square de- 
gree range. Given the realistic response time of an X-ray 
observatory such as IXO, this is the error box relevant 
to X-ray searches for electromagnetic counterparts. 

How easily could such a source be observed? Con- 
tinuously monitoring a whole LISA error ellipse on the 
sky in X-rays does not appear to be a possibility since 
no X-ray telescope is currently planned that has a field- 
of-view of several square degrees and a flux sensitivity of 
10 -15 erg cm -2 s _1 over ~ 10 4 seconds of exposure time. 
A more successful strategy would be to observe a precur- 
sor electromagnetic event, identify a source, and monitor 
it throughout the coalescence. The International X-ray 
Observatory has a planned field of view of 18' (diameter) 
and a flux sensitivity of 10 -15 erg cm -2 s _1 for exposures 
of 10 4 seconds. Thus, approximately 15 frames would be 
required to cover a square degree, which at 10 4 seconds 
for each pointing amounts to approximately two days. At 
the planned sensitivity of the IXO the X-ray outbursts 
could be detected for binaries with orbital periods of a 
few hours and longer, in the frame of the observer. At 
least two subsequent snapshots of the same sub-field of 
view are necessary in order to detect such an outburst. 
The monitoring of the precursor increase in luminosity of 
the disk requires a similar strategy (multiple exposures 
of the same field of view), except the baseline between 
the observations of the same sub-field can be longer than 
in case of the outbursts and at most a few days apart. 
Note, however, that in cases for which the LISA error 
ellipse is smaller than the IXO field of view, the X-ray 
observatory may immediately focus on that part of the 
sky without needing to sweep over a larger area. If a 
candidate with an EM precursor were found, continuous 
monitoring by IXO would be able to detect the flux drop- 
out and subsequent recovery predicted by our models. 

Concurrently, the LISA error box could be monitored 
in the radio band. The Square Kilometer Array is pro- 
jected to have a one square degree field of view and a 
1/iJy sensitivity for a one hour exposure. At 10 GHz, 



1/iJy is 10 -19 erg cm -2 s _1 , so if more than a fraction 
~ 10 -5 of the total emission is in radio waves this will 
be detected easily, as will the change in emission during 
the preceding inspiral. 

If it is possible to identify a counterpart to within a 
few arcseconds, then further followup observations could 
be pursued in other wavebands with yet higher angular 
resolution. We note that if few-arcsecond positions are 
possible then since magnification variations from weak 
gravitational le nsing are dominate d by scales above tens 
of arcseconds ( |Dalal et al.||2003| it might in principle 
be possible to construct detailed shear maps of the re- 
gion and get a moderately refined estimate of the actual 
magnification, although this would require a high surface 
density of background galaxies, in excess of 10 6 deg -2 . 
If this procedure yields better than the us ual ~ 10% pre- 
cision limited by lensing, as discussed by |Jonsson et al. 
(2007), these sources would serve as very precise cosmo- 



logical probes. 



5. CONCLUSIONS 



We have presented a combined analytic and numerical 
analysis of the effects on a circumbinary accretion disk 
of central mass loss caused by prompt gravitational wave 
emission from a binary black hole inspiral and merger. 
Our significant results are: 

1) The primary effect of abrupt mass loss caused by 
the merger of supermassive black holes will be a rapid 
drop in the mass accretion rate across the inner edge of 
the disk. If the mass accretion rate is close to Edding- 
ton and the inner edge can extend without significant 
disruption close to the innermost stable circular orbit in 
the last stages of binary inspiral, then this drop could re- 
sult in large changes in the X-ray spectrum due to, e.g., 
dramatically reduced Comptonization from hot optically 
thin gas inside the ISCO. The flux changes from this ef- 
fect would be detectable with current instruments if the 
source were located precisely enough, but a full survey of 
the likely error ellipses would require an instrument with 
a field of view of at least several tenths of a square de- 
gree and sensitivity in a ~ 10 4 second exposure to fluxes 
below ~ 10 -15 erg cm -2 s _1 . 

2) The expected enhancement of disk emission follow- 
ing a binary merger due to the circularization of elliptical 
orbits is not likely to be observed. While shocks are iden- 
tifiable in our simulations, they contribute little to the 
bremsstrahlung emissivity of the disk. Moreover, such 
shocks are separated by rarefactions that act to reduce 
the brightness of the disk. 

3) In addition to serving as probes of cosmology and 
general relativity, observations of EM signatures from 
merging binary black holes have the potential to tell us 
much about the inner structure of accretion disks. If 
the observed emission from such systems is tied to the 
accretion rate in particular, the duration of the observed 
flux deficit can be used to estimate the viscous timescale 
of the disk. 
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